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Abstract 

In this paper we consider general nearly integrable analytic Hamiltonian systems of one 
and a half degrees of freedom which are a trigonometric polynomial in the angular state 
variable. In the resonances of these systems generically appear hyperbolic periodic orbits. 
We study the possible transversal intersections of their invariant manifolds, which is expo- 
nentially small, and we give an asymptotic formula for the measure of the splitting. We see 
that its asymptotic first order is of the form Ke^e~ a / £ and we identify the constants K, j3, a 
in terms of the system features. We compare our results with the classical Melnikov Theory 
and we show that, tipically, in the resonances of nearly integrable systems Melnikov Theory 
fails to predict correctly the constants K and (i involved in the formula. 

1 Introduction 

In this paper we consider nearly integrable analytic Hamiltonian systems of one and a half 
degrees of freedom. These Hamiltonians are of the form 

h(x,I,T;5) = hoi^ + Sh^xJ^S), (1) 

where 5 <C 1 is a small parameter, (x, r) € T 2 , / 6 U C R and ho and h± are analytic functions. 

When 5 = 0, the Hamiltonian system is completely integrable (in the sense of Liouville- 
Arnold) and / is a first integral of the system. That is, for 5 = the phase space is foliated by 
invariant tori 1 = 1* with frequency vector ui(I*) = (diho(I*),l). The dynamics in these tori 
is quasiperiodic or periodic depending whether difiQ(I*) is rational or irrational. The natural 
question concerning nearly integrable Hamiltonian systems, is how the dynamics changes when 
one adds the perturbation. Namely, which tori persist and which break down, and, what new 
kinds of dynamics appear in the regions where the unperturbed invariant tori break down. 

The persistence of most of the invariant tori (in the measure sense) was shown by Kolmogorov- 
Arnol'd-Moser Theory (see [Kol541 IArn631 lMos 62] ) . The purpose of this paper is to study some 
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particular dynamics in the complementary of this region. Namely, to study the dynamics in the 
region in which the tori that existed for the unperturbed system break down. These regions are 
usually called resonances and they appear when the frequency vector co(I*) of the unperturbed 
system is rationally dependent. This happens for 1 = 1* such that djho(I*) G Q. We call 
n/m : 1-resonance to the region of the phase space I ~ I* such that <9//io(/*) = n/m, n € 7L 
and m G N. Under certain non-degeneracy conditions, in the resonances appear, among other 
invariant objects, low dimensional hyperbolic invariant tori which, for systems of one and a half 
degrees of freedom, are just hyperbolic periodic orbits. These objects have stable and unstable 
invariant manifolds, often called whiskers. The main goal of this paper is to study the possible 
transversal intersection of these invariant manifolds. These transversal intersections ensure the 
existence of a horseshoe, and therefore of chaotic dynamics |Sma65j . 

Even if there is a standard theory to analyze the intersection of stable and unstable invariant 
manifolds, the so called Melnikov Theory (see [Mel63j and [GH83| for a more modern exposition), 
as we will see, it does not apply to this kind of systems. The reason is that we are in a 
singular perturbative setting (see Section |2.1|) . Indeed, as A. Neisthadt proved in |Nei84] . in 
this setting, when the system is analytic, the invariant manifolds are exponentially close to 
each other with respect to the parameter 5. This implies that there is no hope that classical 
perturbative techniques, such as Melnikov Theory, will allow to study the possible intersection 
of these manifolds. 

The problem of exponentially small splitting of separatrices was first pointed out by Poincare 
[Poi90l §19] and has drawn considerable attention in the past decades (see [HMS88, ISMH911 
IDS921 IFon931 IGG941 iGeiM IFon951 ISau951 IDG.7S971 IDS971 IGel97al IGel97bl lTFe97l IGGM991 
iGeinni ISauOll IOSS031 IDGS041 IBFMI IBFOKI iBalMl lUIiMl IGOS101 IBFGS111 IGSTT] ). Neverthe- 
less, most of these results cannot be applied to the resonances of nearly integrable systems since 
they assume an artificial smallness condition for certain terms in the perturbation hi. The only 
ones that can be applied to this kind of systems are |Tre971 IGelOOl IOSS031 IOh061 IGOSlOj , but 
they only deal with particular examples. The present paper is strongly related to |BFGS11| . 
In that paper, the authors study the exponentially small splitting problem for slightly differ- 
ent systems (see the paper for the exact hypotheses). Nevertheless, one can easily see that 
Hamiltonians (pQ) with ho (I) = I 2 /2 and h± independent of I studied in the : 1 resonance fit 
the framework of [BFGSli] and therefore, the results obtained in that paper can be applied 
to these systems. In the present paper, using the tools developed in jBFGSll] , we study the 
exponentially small splitting problem for more general Hamiltonians of the form ([1]) close to any 
resonance (in particular allowing general ho and general dependence on / of hi, see Section [2.21 
for the exact hypotheses). Nevertheless, as happened in [BFGS1T] . we have to assume that h\ 
is a trigonometric polynomial in x. This hypothesis has been assumed in all the previous works 
since as far as the author knows, there are not known tools to deal with more general (analytic 
in x) perturbations. 

The structure of the paper goes as follows. First, in Section [2] we study the splitting of 
separatrices for Hamiltonians of the form ([1]) in the : 1 resonance. In this section we state the 
hypotheses we need to impose to Hamiltonian ([TJ (Section I2.2p and we state the main results 
(Section I2.3p . In Section l2~4l we apply the obtained result to a particular (and paradigmatic) 
example. From the results obtained for the : 1 resonance, in Section [3] we deduce analogous 
results for all the other n/m : 1-resonances. Finally in Section H] we prove the results stated in 
the previous sections. 



2 



2 Results for the : 1 resonance 



2.1 The rescaled Hamiltonian System in the resonance 

To facilitate the exposition, in this section we focus on the : 1-resonance. Namely we consider a 
region of the phase space I ~ I* such that I* satisfies that the corresponding frequency vector is 
lo(I*) = (djh()(I*), 1) = (0, 1). Any other resonance can be reduced to this one after a change of 
variables as will be seen in Section[3j After a translation in /, one can assume that co(I*) = (0, 1) 
takes place at I* = 0, which implies that <9//io(0) = 0. Then ho is of the form 

h (I) = /i 02 y + 9(1) with g(I) = O (/ 3 ) . (2) 

We assume that /102 ^ 0, which guarantees that the resonance is non-degenerate (see Section 
12.21 for the exact hypotheses we need). Then, without loss of generality, we can also assume 

h 02 > 0. (3) 

The dynamics of the perturbed system around this resonance, is better understood if one 
performs the rescaling 

—y and r = J_ , (4) 
no2 Vh-020 

since then we magnify the size of the resonant zone, so that its width becomes of order one with 
respect to 5. To simplify the notation, we take 



(5) 



as a new parameter. Then, one obtains a Hamiltonian System with Hamiltonian 

,2 



where 



H (x, y, t -;e^j=^- + ^G(ey) + V(x) + F (x,-\ + R (x, ey, l;eY (6) 



G(I) = h 02 g ( j- ) (7) 
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R(x, I, r; e) = h 1 [ x,-^,t;^— ) - h^x, 0, r; 0). 



V(x) = (hi{x, 0, r; 0)) = — / h(x, 0, r; 0) dr (8) 

27r Jo 

F{x,t) = /n(rE,0,T;0) - (h^x, 0, r; 0)) (9) 
ho2 ' ' h.02 , 

Note that there are terms of different size. The terms y 2 /2 and V(x) are independent of the new 
parameter e and have become the main terms of the rescaled Hamiltonian. The term F(x,t/e) 
is also of order one with respect to e but it is rapidly oscillating and has zero average. Thus, as 
it is well known from Averaging Theory (see for instance |AKN88j ). this term is in fact causing 
effects of order 0(e) to the system. More precisely, one can perform an e-close to the identity 
symplectic change of variables so that this term becomes of order 0(e). Finally, the terms 
e~ 2 G(ey) and R(x,ey,t/e;e) are of order 0(e). 
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Taking into account these different sizes, one can split Hamiltonian ([6j) as a new unperturbed 
Hamiltonian plus a new perturbation. Namely, one can rewrite it as follows. 



H (x,y, ^\ = H (x,y) + /j,Hi fx,y, pe^j (11) 

with 

H (x,y) = £ + V(x) (12) 

Hi(x,y,r;e) = F(x,r) + \c(ey) + R(x,ey,r;e). (13) 

Note that here [i is a fake parameter since we are interested in fx = 1. Nevertheless, we write 
it to clearly stress that now we have a new unperturbed system Hq (that is, H with \i = 0) 
and a perturbation \iH\. Moreover, stating the results including the parameter \x makes them 
simpler to compare them with the previous ones which usually also include this parameter. 
Indeed, the first results in the area assumed some smallness condition on this parameter \x (see 
[HMS881 IDS921 IGel941 IDS971 IGel97bl lBF04l [BF05]). In the case fj, = 1 there are fewer results. 
In |Tre97l IGelOOj IGUSlOj the authors study particular examples, which in particular satisfy 
G = and R = 0. In [Bal06l IBFGSli] the authors study systems of a slightly different form. 
Nevertheless, one can see that general Hamiltonians of the form (jlip with G = and R = fit 
their framework. Here we generalize these results to Hamiltonians (|11|) with general G and R. 
In particular we show that the functions G and R contribute to the main order of the difference 
between the invariant manifolds and, in particular, if one omits them one obtains a wrong size 
of the asymptotic order. 

The equations associated to the Hamiltonian (jlip are 

x = y + fidyHi (x,y, z ;e 



y 



-V'(x) - nd x H x fx, y, p e'j 



(14) 



From now on, we call unperturbed system to the system defined by the Hamiltonian Hq and we 
refer to Hi as the perturbation. 



2.2 Hypotheses 

We devote this section to state the hypotheses needed for the main results concerning the 
resonance 0:1. Note that all the hypotheses refer to the integrable Hamiltonian ho and the 
first order in 5 of the perturbation hi in (pQ). Namely, the 0{5 2 ) terms in (pQ) do not play any 
role. Some of them are stated using the rescaled hamiltonian pip . 

HP1 There exists r > independent of 5, such that the Hamiltonian h in ([1]) is analytic for 
(x,I,t) E T x B(r) x T, where B{r) is the ball B(r) = {I € C : |/| < r}. Moreover, h 
satisfies 

d]h (0) > 0. 

HP2 The original Hamiltonian is of the form ([1]) and h\ is a trigonometric polynomial of degree 
M > 0. 
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HP3 The potential V in (JSj) is a trigonometric polynomial of the same degree as h\ in ([I]). 
Namely, V has degree M. 

HP4 The system associated to Hq has a hyperbolic critical point at (0, 0) with eigenvalues 
{A, —A} with A > 0. Equivalently, the potential V defined in ([8]) satisfies 

A 2 

V(x) = — —x 2 + O (x 3 ) as x ->■ 0. 

HP5 The stable and unstable invariant manifolds of the critical point (0, 0) of the system 
associated to the Hamiltonian Hq coincide along a separatrix. 



We denote by (qo(u),po(u)) a real- analytic time parametrization of the separatrix given by 
Hypothesis HP5 with some chosen (fixed) initial condition. It is well known (see |Fon95j ) 
that there exists a > such that the parametrization (qo(u) , po(u)) is analytic in the complex 
strip {|Imn| < a} and not in any wider strip. Next hypothesis deals with the behavior of this 
parameterization at the boundary of this strip. 

HP6 There exists a real-analytic time parametrization of the separatrix (qo(u) , po(u)) analytic 
on {|Imu| < a} such that the only singularities of (qo( u )iPo( u )) i n the lines {Imu = ±a} 
are at u = Ha. 

Note that the strong assumption in Hypothesis HP6 is having only one singularity at each 
boundary. Once there exists only one, one can just make a time translation to locate it at the 
imaginary axis. 

Together Hypothesis HP2 and HP6 imply that qo(u) has logarithmic singularities at ±ia 
of the form qo(u) ~ \n(u =p id) (where we take different branches of the logarithm whether we 
are close to +ia or — ia: we take arg(u — ia) € (— 3ir/2, vr/2) and arg(u + ia) G (— ir/2, 3ir/2) 
respectively). In this case, one can see that, if u £ C, \u =F ia\ < u, 

cos(q (u)) = °± (l + O ((« =F iafl M 

*"*'-»= (^JM i+o (^» )2/M )) (15) 

'»w = (^( 1+0 ((»^ i °) 2/M 

with arg(u — ia) £ (— 3tt/2, n/2) and arg(u + ia) £ (— tt/2, 3tt/2) if we are dealing with the 
singularity +ia or —ia respectively (see Section 2.1.3 of |BFGSli] ). We also have that 

C+ = CZ = ±i^- 

Remark 2.1. Hypothesis HP1 can be weakened to 9f/io(0) ^ and both the positive and 
negative case can be proved analogously. To simplify the exposition we assume djho(0) > 
(the other case can be deduced from it just reversing time). Note that this condition has already 
been stated in (|3j) and has been used for the rescaling and obtention of the model close to the 
resonance ([T]). 
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Remark 2.2. The critical point (0,0) has invariant manifolds thanks to Hypothesis HP4- 
Therefore, the actual assumption in Hypothesis HP 5 is that they coincide along a separatrix. It 
is easy to see that the energy levels of Hq are compact sets and therefore, the only possibilities 
is that the invariant manifolds of (0, 0) coincide with either the invariant manifolds of the same 
or another critical point. Therefore, HP5 rules out the case of heteroclinic connections. The 
existence of heteroclinic connections is a non-generic phenomena, since after a small perturba- 
tion one can make that all the critical points of Hq belong to different energy levels and thus 
impossible to be connected by heteroclinic connections (see Section \2.3.3\) . 

Remark 2.3. The separatrices, whose existence is guaranteed by Hypothesis HP5 can be of two 
different topological type. Indeed, note that the hyperbolic critical points of the system associated 
to Hamiltonian (|12p are of the form (x*,0), where x* are non- degenerate maxima ofV. Then, 
if the maximum is global, the associated separatrices are a graph over the variable x, as happens 
for the pendulum, which has only one maximum (see Figure^). If the maximum is local but 
not global, the associated separatrices form a figure eight as can be seen in Figure UJ Therefore, 
if V has more than one maximum there always coexist graph-like and figure eight separatrices 
associated to different critical points. The results obtained in this paper can be applied to both 
types of separatrices as long as Hypotheses HP1-HP6 are satisfied. 



V 




Figure 1: Phase portrait of the system associated to Hamiltonian (|12p when the potential V has 
only one local (and thus global) maximum, which after a translation in x can be placed at the 
origin. Note that this maximum has to be non-degenerate to ensure that the critical point is 
hyperbolic. See Figure [2] for the phase portrait when V has more local maxima. 

2.3 Main results for the : 1 resonance 

We devote this section to state the main results about the : 1 resonance. First we state 
Theorems 12.41 and 12. 5| which give the main results referred to the rescaled Hamiltonian (|lip and 
then, from them, we deduce Theorem 12. 6| which states the main results for the resonance : 1 
in the original variables. That is, we give it in terms of Hamiltonian ([TJ. In Section Owe deduce 
from Theorems 12.41 and 12.51 analogous results for general resonances n/m : 1. 

By Hypothesis HP4, the system associated to the Hamiltonian Hq in (|12|) has a hyperbolic 
critical point at the origin. Next theorem ensures that the hyperbolic critical point of the 



G 



y 




Figure 2: Phase portrait of the system associated to Hamiltonian (|12p when the potential V 
has one global and one local maximum. The global one corresponds to the critical point whose 
separatrices are a graph over x. The local (and non-global) critical point is the one having figure 
eight separatrices. If V has more local maxima, one obtains more nested figure eight separatrices 
in the phase space 

unperturbed system becomes a hyperbolic periodic orbit for the full system (|14|) . which is close 
to the origin. 

Theorem 2.4. Let us assume Hypotheses HP1 and HP4- Then, given any fixed value [i > 0, 
there exists £q > such that if e G (0, £q), the system associated to (|14p has a hyperbolic periodic 
orbit (x p (t / e) , y p (t / e)) which satisfies that, for r G M, 

|x p (r)| <^Me 2 
\Vp(t)\ < K\n\e, 

for a constant K > independent of e. 

The proof of this theorem follows the same lines as the proof of Theorem 2.2 in [BFGSll] 
(see also [Fon93Hr5n551ID557| ). 

The next step is to study the stable and unstable invariant manifolds of the periodic orbit 
(x p , y p ). In the unperturbed case (that is fx = 0) we know that they coincide along the separatrix 
(qo,Po) given in HP6. We will see that when \i ^ they tipically split and their transversality 
is exponentially small. 

To measure the splitting of the invariant manifolds we consider the 27re-Poincare map Pt in 
a transversal section T, to = {(x,y,to); (x,y) G M 2 }. This Poincare map has a hyperbolic fixed 
point (x p (to / e) , y p (to / e)) , whose stable and unstable invariant manifolds are curves. 

As Pt is an area preserving map, we measure the splitting giving an asymptotic formula 
for the area of the lobes generated by these curves between two transversal homoclinic points. 
Moreover, by the area preserving character of Pt , the area A of these lobes does not depend on 
the choice of the homoclinic points. Other quantities measuring the splitting, as the distance 
along a transversal section to the unperturbed separatrix, or the angle between these curves at 
an homoclinic point, can be easily derived from our work. 

Theorem 2.5. Let us assume Hypotheses HP1-HP6. Then, given any fixed \x, there exists 
Sq > such that if £ £ (0,Eq), the area of the lobes between the invariant manifolds of the 
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periodic orbit given in Theorem \2.J\ is given by 



A a\ 1—1 — — +itlm 6 In - 

A = 4\fi\e e e p e 



iC(tM) 



+ o 



hid 



(16) 



where f(/i) is an entire analytic function, C(fi) is an entire analytic function such that C(fj.) 
O(n) and b £ C is a constant defined as 



b = 2g 3 C^ 



(17) 



where C + has been defined in (fTS"]) and g 3 is the degree 3 Taylor coefficient of the function G in 
jZD, namely, G{I) = g 3 I 3 + 0(I 4 ). 

The proof of this theorem is given in Section^ Recall that we were primarily interested in 
the case fx = 1, in this case the asymptotic formula for the area becomes 

1 



A = Ae~ 1 ~ lmb e-e 



iC(l) 



+ o 



hid 



(18) 



Theorem 12.51 can be rephrased to give an analogous result for the system associated to the 
original Hamiltonian h in ([I]). 

Theorem 2.6. Let us assume Hypotheses HP1-HP6. Then, there exists 5$ > such that if 
5 £ (0, So), the system associated to Hamiltonian has a hyperbolic periodic orbit (x p (t), I p (t)) 
which satisfies 

Mr) | < K5 
\I p (r)\ < K5. 

Moreover, the area of the lobes between its invariant manifolds is given by 



A = A(h 



Im b Im b 



02, 



1 



ln<5| 



(19) 



where /102 = d]ho(0) (which is positive by Hypothesis HP1), = f(l)e tC ^ and b, /(l) and 
C(l) are the constants given by Theorem \2. 5[ 

Note that to deduce this theorem from Theorem 12.51 it is enough to undo the changes ([!]) 
and (0). 



2.3.1 Some remarks 

In the first order of formulas (|18p and (|19p there appear three constants /(l), C(l) and b. We 
devote this section to make some comments about them. 

• The constants /(l) and C(l) play a similar role. Nevertheless, the functions /(//) and 
C(n) have a significantly different origin, which can be seen in Section [U Let us just say 
here that C(n) can be obtained by classical perturbative methods, in the sense that it can 
be computed in terms of the first orders in e of the funtions R(x, ey, r; e) and e~ 2 G(ey). 
Instead, f(n) has a completely different origin. It comes from the study of the inner 
equation associated to system (fT4|) . which is an equation independent of e and was studied 
in |Bal06j (see Section [J] for its definition) . In Section U we will see that it depends on the 
full jet in / of the funtions G{I) and R(x, I, r; e). Namely, it depends on all the orders in 
e of R(x, ey, r; e) and e~ 2 G(ey). Note that this implies that any finite order truncation in 
e of the Hamiltonian H does not predict correctly the area of the lobes. 
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• To ensure that the area of the lobe is indeed given asymptotically by formula (|18p . one 
has to check that /(l) 7^ 0. With the techniques used in this work, as far as the authors 
know, there is not any analytical way to check this condition. Nevertheless, one can use 
the mentioned inner equation to verify this condition numerically |SV09| . 

• The constant b affects the polynomial term in the formulas for the area of the lobes (|18p 
and (|19p . Generically satisfies b 7^ 0, since generically <?3 7^ 0. Nevertheless, it vanishes 
when 53 = 0. In previous works (see |Tre97l IOli06] ) . the models considered satisfied G = 
and thus this coefficient did not appear in the formulas provided by those works. This 
term was first detected in |Bal06l IBFGSll] . Note that, as it happens for the function 
C(fj,), it can be obtained by the first order in G, that is, one does not need the full jet of 
this function to compute it, as happens for /(//). 

2.3.2 Comparison with Melnikov Theory 

If one considers the system associated to Hamiltonian (jlip , one can apply the classical Melnikov 
Theory, which is just to make a perturbative approach taking \x as a small parameter. Recall 
that we are interested in the case \i = 1 and not in arbitrarily small /i. Nevertheless, this was 
the first approach to study the phenomenon of exponentially small splitting of separatrices since 
turns out to be the simplest case. Indeed, using this method and taking [i small enough, one 
can see that the area of the lobes is given asymptotically in \i as 

A = Ae- l e-~e (\f \ ^ + O (f, 2 )) , (20) 

for certain constant /o € C. 

Thus, if one wants to assure the validity of the first order given by this formula, it has to be 
assumed that fi is exponentially small with respect to e. Nevertheless, it is a well known fact 
that Melnikov gives the correct first order for a wider range in \x (see [DS971 IBFGSll] ) . Using 
Theorem 12.51 one can see for which range of \x Melnikov indeed gives the correct first order and 
for which range it does not. In [Bal06j . it is seen that the function f(n) appearing in Theorem 
12.51 and the constant /o involved in the Melnikov prediction (|20p are related by 

m = f +o(fj). 

Thus, comparing formulas (|16|) and (|2U|) . one can see that in the general case b 7^ 0, the Melnikov 
prediction is correct provided 

1 



In the non-generic case b = 0, which takes place for instance when G = 0, it is enough to require 
fj,< 1. 

In the resonances of nearly integrable Hamiltonian Systems, the parameter fi satisfies fi = 1 
and therefore Melnikov does not predict correctly the area of the lobes. Indeed, comparing 
formulas (|16p and (|20p . one can conclude that Melnikov only predicts correctly the coefficient 
in the exponential, but fails to predict correctly both the polynomial term and the constant 
coefficient. 
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2.3.3 On the genericity of the hypotheses and the results 

We devote this section to make some remarks on the genericity of Hypotheses HP1-HP6 stated 
in Section [2.21 assumed in Theorem 12.51 and 12.61 an d about the genericity of the results stated 
in these theorems. 

Let us consider a constant r > 0. Then, it is clear that HP1, which is the only condition 
needed for the integrable system ho in ([TJ, is satisfied for a generic ho in 

X = {/ : {\I\ <r}cC4C: real-analytic, WfW^ < oo} , (21) 

where || • ||oo denotes the classical supremmum norm. 

Concerning the hypotheses imposed on the perturbation h\ in ([TJ), let us fix also constants 
> and a > 0. Then, one would like to obtain their genericity in the space 

y = {/ : T a x {|/| < r} x T CT x {|e| < e } -> C : real-analytic, H/IU < oo} , 

where T CT = {x G C/Z : |Imx| < a}. Nevertheless, functions hi £ y satisfying Hypothesis HP2 
are clearly non generic since they are trigonometric polynomials in x. However, one can ask 
about the genericity in the smaller space of perturbations which are a trigonometric polynomial 
in x of degree M > or less. Thus, we consider the space 

( M 

y M = \f= f [k] (I,r;e)e ikx : x {|/| < r}xT a x {|e| < e } -> C; 

k=-M (22) 

real-analytic, ||/||oo < oo 

which is a Banach space contained in y. By definition, HP2 is generic in this space. Taking 
into account the definition of V in terms of hi in ([8]), one can easily see that HP3 is also generic 
in 3^m- The same happens for HP4 since, as explained in Remark 12. 3( the hyperbolic critical 
points of the system Hq in (|12p correspond to non-degenerate maxima of the potential V and 
it is a well known fact that for a generic V all its extrema are non-degenerate. Moreover, in 
Remark 12. H we have already explained that having a homoclinic connection for a hyperbolic 
fixed point is also generic. Thus we can conclude that Hypotheses HP2-HP5 are generic. 

It only remains to check HP6. It is clear that this hypothesis is open in J^m- Let us check 
its genericity. As we have seen in Remark 12.31 there are two types of separatrices: the ones that 
are graphs with respect to x and the ones that form a figure eight (see Figures [Tj and [2|) . Let us 
fix a concrete separatrix and check the genericity of Hypothesis HP6. 

In the figure eight case, the separatrices are reversible with respect to the involution $(x, y) = 
(x, —y). Then, Hypothesis HP6 is open but not generic. Indeed, if for a concrete potential, a 
figure eight separatrix has singularities at points u = ±a ± ai with this structure cannot 

be broken by any small perturbation of the potential due to the reversibility. 

Namely, for this kind of separatrices in 3^w" there is an open set where Hypothesis HP6 is 
satisfied and a different open set where the separatrix has two singularities at each connected 
component of its strip of definition. 

For the graph separatrices, as a far as the author knows, there are no results in the litera- 
ture about the genericity of Hypothesis HP6. Nevertheless, since in this case the separatrices 
generically do not have any symmetry, one would expect that if a separatrix has singularities 
at u = ±a ± ai, one can modify slightly the potential V so that the singularities have different 
imaginary part. This fact, prompt us to state the following conjecture. 
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Conjecture 2.7. Let us consider a Hamiltonian of the form 

H (x,y) = ^ + V(x). 

Then, for a generic trigonometric polynomial V , the associated Hamiltonian system has a hy- 
perbolic critical point, such that its invariant manifolds coincide forming a separatrix which is 
a graph with respect to x and satisfies Hypothesis HP6. 

Recall that for a generic potential there always exists a hyperbolic critical point having 
graph separatrices, since they correspond to a non-degenerate global maximum of V . Thus, if 
this conjecture is true, the assumed hypotheses are satisfied by a generic ho G X (see ([2"T]) ) and 
a generic h\ € 3^m (see (|22p ). Nevertheless, having generic hypotheses on the Hamiltonians ho 
and h\ does not imply directly that we can deduce from our results that generically the invariant 
manifolds split transversally and that the area of the associated lobes is given by the asymptotic 
formula (|19|) with non-vanishing first order. Indeed, to obtain this result one would need to 
prove that /(l) 7^ for a generic set of Hamiltonians. 

Let us show that assuming Conjecture 12.71 and adding an extra generic hypothesis this is 
the case (under a suitable definition of genericity). Indeed, as stated in Theorem 12.51 the 
function f(fi) is analytic. Moreover, recall that /(/i) = fo + C(/u) where fo is the constant 
given by Melnikov Theory (see Section T2.3.2p . This constant fo comes from certain residue that 
generically does not vanish, which implies that generically /o / (see [DS971 IBF04j for a more 
detailed explanation). Then, since / is analytic, under this extra hypothesis / only vanishes 
for a discrete set of values of /i and thus for a generic set of Hamiltonians and a generic value 
of /i, we have that /(/i) 7^ 0. Now we show that this fact is sufficient to show that generically 
f{l) 7^ 0. The openness of this condition is clear. Therefore, it is enough to check that it 
is dense. Let us consider ho and h\ satisfying all the stated hypotheses (and that fo satisfies 
fo 7^ 0) such that /(l) vanishes. Since f(/j,) only vanishes for a finite set of points, there exists 
fx* arbitrarily close to 1 such that /(/U*) 7^ 0. Notice that modifying slightly \i means modifying 
slightly Hi, which is equivalent to modifying slightly ho and hi in (TjQ). Thus, for generic pair 
(ho, hi) G X x ym (with respect to the product topology), we have that /(l) 7^ 0. 

Therefore, we can conclude that, if we fix M > and assume Conjecture 12.71 then for a 
generic pair (ho, hi) £ X x 3^m of Hamiltonians, there exists £0 > such that for any e S (0, £q) 
the Hamiltonian (pQ) has a hyperbolic periodic orbit at the resonance whose invariant maniolds 
are graphs, intersect transversally and the area of the lobes is given by the asymptotic formula 
(|19p of Theorem 12. 6| whose first order does not vanish. 

Finally, let us point out that, if one takes into account the rescaled formula (|16p . we can 
ensure that for a generic set of Hamiltonians (in the sense just explained), Melnikov Theory 
does the wrong prediction for the area of the lobes both in the polynomial and constant terms 
(see Section [2. 3. 2p . 

2.3.4 The parabolic case 

The generic Hypothesis HP4 assumes the that the Hamiltonian Ho in (|12p has a hyperbolic 
critical point at the origin. Nevertheless, the results obtained in Theorem [23] can be generalized, 
under suitable hypotheses, to Hamiltonian Systems which have a parabolic point. Indeed, let us 
assume that Ho has a parabolic point at the origin which has invariant manifolds that coincide 
along a separatrix, as happens, for instance, for 

H (x,y) = y - (cosx - l) 2 . 
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Having a parabolic critical point at the origin, is equivalent to ask that the potential V(x) (see 
(J5D) satisfies 

V(x) = v k x k + O (x k+1 ^ as x -> 

for certain k > 3 and 6 R. Then, to avoid bifurcations, the extra hypotheses that one has 
to assume are to guarantee that, when we add the perturbation, the critical point preserves 
its parabolic character and that its invariant manifolds are preserved. To this end one has to 
impose the following two conditions. The first one is that the function F in (J9J) satisfies 

F(x, r) = a(T)x m + O (x m+1 ) as x -»• 

with 2m — 2 > k. The second condition is the analogous one for R in (jlOp . Indeed, if we 
Taylor-expand R with respect to / around I = we get 

R(x, I, r; e) = R x (x, r; e)I + O (I 2 ) 

Then, the second additional hypothesis is that R± satisfies 

Ri(x, r; e) = /3(r; e)x n + O (x n+1 ) as x -> 

with n > fe/2. Under these two extra hypotheses, Theorem 12.51 is also valid for parabolic points. 
The proof of this fact was done for slightly different Hamiltonian Systems in [BFGSll] (see also 
|BF04j ). and it can be easily adapted to our setting. Note that these two extra hypotheses are 
satisfied in the natural case that F and R are of the same order as V, that is if F = 0{x k ) and 
R = 0(x k ). However, the hypotheses we need to require are considerably weaker. 



2.4 An example 

In this section we apply Theorems 12.51 and 12.61 to a particular example. We consider the Hamil- 
tonian 

I 2 

h(x,I,r;S) = — +r]I 3 + 5(1 + al) (1 + sinr)(cosx - 1), rj,aER, (23) 

and we study its splitting of separatrices for the resonance 0:1, which takes place at / = 0. We 
have chosen this model for two reasons. Firstly, (most of) the constants involved in the splitting 
formula can be computed explicitly and thus we can analyze it in great detail. Secondly, it shows 
all the modifications that the splitting formula of Theorem 12.51 presents in comparison to the 
Melnikov prediction. Moreover, it can be compared with the results in |Tre97j . which considers 
the same model with ij = and a = 0. 

Note that for this model /102 = 1, and thus the rescaling is simply y = \/SI and e = y/S. 
Then, it is easy to see that the rescaled model is given by 

H (x,y,^ = H (x,y) + H! (x,y,?j (24) 

with 

y 2 

H (x,y) = — + cosx - 1 
Hi(x, y, t; e) = sinr(cos x — 1) + aey(l + sin r) (cos x — 1) + rjey 3 
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Thus, the functions V, G and F and R associated to this system (see ((HJ), ©, ([9|) and (fTU|) ) are 
given by 

V(x) = cos x — 1 
G(I) =r ? / 3 

F(x,t) = sin r (cos x — 1) 
R(I,x,t) = al(l + sinr)(cosa; — 1) 

The Hamiltonian Hq is the Hamiltonian of the classical pendulum and therefore Hamiltonian (|23|) 
satisfies Hypothesis HP4 and HP5. The upper separatrix of the pendulum can be parameterized 
as 

j(u) = I 4arctan (e u ) , — 

\ cosh u 

and its singularities are located at u = iir/2 + k?ri, k € Z. Therefore Hypothesis HP6 is also 
satisfied. It is easy to see that Hamiltonian (|23j) also satisfies Hypotheses HP1-HP3 and thus 
one can apply Theorem 12 .51 Note also that the periodic orbit at (x,y) = (0, 0) is preserved when 
one adds the perturbation and, therefore, Theorem 12.41 is not needed. 

We have that g$ = rj and C+ = — 2i and therefore, using formula (|17|) . the constant b is given 
by b = —Arji. Then, the asymptotic formula for the splitting for system (|24j) reads 



A = 4e- l+Ar) e~Te 



iC(l) 



As we have already explained, the constant /(l) cannot be computed with our methods and 
therefore, we cannot ensure whether it vanishes or not. Even if it is not written explicitly, both 
/(l) and C(l) depend analytically on the parameters r] and a. 

If we take rj = and a = we recover from (|25p the formula given in |Tre97] . As we have 
said only when 77 = 0, the constant b vanishes and then, the degree of the polynomial term 
coincides with the Melnikov prediction. Note also that one could add higher orders in / both 
in ho and h\. Then, in the splitting formula (|25j) only would change /(l) and C(l) but b would 
remain the same (see Q17])). 

Undoing the rescalings, one can deduce from (|25|) the formula for the original system (|23p . 
which reads 

k« eiC(I) l +0 (n^ 



A = 45 2rt e 2V5 



3 Results for a general resonance 

The results stated in Theorem 12.61 are focused on the resonance with frequency vector uj = (0, 1). 
Nevertheless, as we will show in this section, from Theorem 12.51 one can deduce results for any 
other resonance. Indeed, we will reduce any resonance to the resonance 0:1. 

As it is well known, a resonance takes place at / = I* if diho(I*) G Q. Let us assume, 
without loss of generality that the resonance takes place at I* = and then <9//io(0) = n/m 
where n and m are coprime integers. We consider a Hamiltonian ho + Shi of the form (pQ) and 
we show that there exists a change of variables such that the new Hamiltonian ho + Shi satisfies 
djho(0) = 0. Namely, we transform the resonance m/n : 1 into 0:1. Indeed, this change is 
simply given by 

~ n 

x = x r, (26) 

m 
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and the transformed Hamiltonian reads 

h(x, I, t ; 5) = ho(I) + Shi (x, I, t ; 8) 

with 

— . 77 

h (I) = ho(I) - -I 
m 

hi (x, I, r; 5) = hi (x r, /, r; 5] . 

\ m J 

Note that with this change the perturbation becomes 2-7rm-periodic in r. Then, one can apply 
the rescaling in @. Nevertheless, to proceed as for the : 1-resonance, we want that after the 
rescaling the perturbation H\(x, y, t/e; e) is 2tte periodic in time (in the variable t). To this end, 
we just need to define as a new parameter 



m 



instead of ©. 

Note that ho2 = ^02 since the change (|26p only modifies the first coefficient of the Taylor 
series of ho. 

After the rescaling we get a system of the form (jlip with the modified functions V, F, G 
and R defined as 



V(x) 



1 

2^ 



2tt 



hi ( x H r, 0, r ) dr 

m 



F (x, t) = hi ( x H r, 0, r 

m 



V(x) 



i? (x, r; e) = /ti ( x -\ r 



n 



I e 2 
,r: 



m ' /io2 r ^' ' ho2m 2 
I 



hi [ x H r, 0, r; 



777 /i, 



02 



Then, from Theorem [23] we can deduce an asymptotic formula for the splitting of separatrices 
for any resonance. Nevertheless, note that now the parameter e is given by e = ra^Jh^S and 
therefore the necessary smallness condition on the perturbation, e G (0, £0), depends on m, 
namely on the resonance: the higher the resonance is, the smaller 5 is needed. 

Theorem 3.1. Let us assume Hypotheses HP1-HP6 and let us fix a rational number n/m, 
where n andm are coprime integers. Then, there exists 5$ > such that if 5 G (0, 8q), the system 
associated to Hamiltonian ([!]) has a hyperbolic periodic orbit (x p (t), I p (t)) which satisfies 



n 



x p( t ) T 

777 



< Km 2 5 



\I p (r)\ < Kb. 

Moreover, the area of the lobes between the invariant manifolds of this hyperbolic orbit is given 
by 



hub 



A = 4h 



Im b 



02 



S 2 e m V^o2<5 



:c(i) 



+ 



1 



ln(m5V2)| 



where b, /(l), C(l) are the constants given by Theorem \2. 5\ and ho2 = d]ho(0), which is positive 
by Hypothesis HP1. 
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4 Description of the proof of Theorem 12.5 



The proof of Theorem 12.51 follows the same lines as the proof of Theorem 2.7 in |BFGS11] . In 
this section, we describe the proof focusing on the main differences from [BFGSli] . 

4.1 Basic notations 

First, we introduce some basic notations which is used in this section. We denote by T = M./ (2irZ) 
the real 1-dimensional torus and by 

T a = {t G C/(2vrZ); |Imr| < a} , 

with a > 0, the torus with a complex strip. 

Given a function h : D x — > C, where D C C is an open set, we denote its Fourier series 

by 

h(u,r) =Y,h [k \u)e ik 

and its average by 

(h) (u) = hJ® (u) = — [* h{u, r) dr. 



Akr 







We fix also some notation for the functions F, G and R defined in ([9]), and ()10p . Thanks 
to Hypothesis HP2 the functions F and R are trigonometric polynomials of degree at most M. 
Then F, which by construction has zero average, can be written as 



M 

ikx 



F(x,t) = ^2 a k(T~)e 



k=-M 



for certain analytic periodic functions : T — > C with zero average. We write the function R 
in (dHD as 

R(x,ey,T;e) = ^ c k i(r;e)e ikx (ey) 1 



k=-M...M 
1>1 



for certain analytic functions cti ■ T — > C. Finally, we write G in (J7]) as 



G{ey) = ^2g k {ey) k 

k>3 



for certain constants g k G R. 

By Hypothesis HP1, the Hamiltonian H in (jlip is analytic in r = t/e. By the compactness 
of T, actually there exists a constant <jq such that H is continuous in T ao and analytic in T CT0 . 
From now on, we fix < a < <jq. 

Throughout the proof of Theorem 12.51 we will use the analyticity in \i. We fix an arbitrary 
value fio > 0. Even if we do not write it explicitly, all functions which appear in this paper will 
be analytic in ^ S B(hq) = {fi G C : < /j,q}. 

From now on, we work with the fast time r = t/e. Then, denoting ' = d/dr, we have the 
system 

x' =s(y + nd y Hi(x,y,T;£)) , . 

y > =-s(V'(x) + t id x H 1 (x,y,T;e)). [ '> 
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4.2 The periodic orbit 

The next theorem states the existence and useful properties of a hyperbolic periodic orbit close 
to the origin of the perturbed system. 

Theorem 4.1. Assume Hypotheses HP1 and HP4 and fix any hq > 0. Then, there exists Sq > 
such that for any \fj,\ < hq and e € (0, Eq), system ([77]) has a 2ti -periodic orbit (x p (t), V v {t)) : 
To- — > C 2 which is real- analytic and satisfies that for r £ T a , 

\x p {t)\ < b \n\e 2 
\y P {r)\ < b \fi\e, 

where bo > is a constant independent of e and fi. 

The proof of this theorem follows the same lines as the proof of Theorem 4.1 in [BFGS11]. 
Note that this theorem gives the same result as Theorem 12.41 The only difference is that it 
proves the existence of the periodic orbit in an analytic strip around the real axis. This is 
necessary to deal later on with analytic parameterizations of the invariant manifolds. 

Remark 4.2. The periodic orbit (x p (r), y p (r)), the Hamiltonians H and Hi, which will be 
defined in (|29p and (|30|) respectively, and the functions involved in these definitions depend on 
the parameters \x, e. From now on, we omit this dependence explicitly. 

Once we know the existence of the periodic orbit, we perform the time dependent change of 
variables 

q = x- x p (t) 

p = y- y P (T~) 

which transforms system (|27p into a Hamiltonian system with Hamiltonian function eH(q,p, r): 

H(q, P , r)=^ + V(q + x p {r)) - V (x p (t)) - V (x p (r)) q 
+ /j,Hi(q,p,r) 

Hi (q, p, t) = e- 2 G(ep, r) + F(q, r) + R(q, ep, r) (30) 

F(q,r) =F(q + x p (r),T) - F(x p (t),t) - d x F(x p (r),r)q 

R(q,£P,t) =R{q + x p (r),ep + ey p (r),r) - R(x p (t) , ey p (r) , t) 
- d x R(x p (T),ey p (T),T)q - d y R(x p (T),ey p (T),T)ep 
G(ep, t) =G {ep + ey p (r)) - G (ey p (r)) - eG' {ey p (r)) P 

We have added the terms V (x p (t)), F(x p (t),t), R(x p (t) , sy p (r) , r) and G {ey p (r)) for conve- 
nience. Note that they do not generate any term in the differential equations associated to H. 
Since |x p (r)| = O (fJ>£ 2 ), we can split the function F as 

F(q,T)=F 1 (q,T)+F 2 {q,T), 



where 



with 
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where 

F 1 (q,r)= Yl a k {T)(e ki i-l-kiq) 
k=-M...M 

and F2(q,r) is the remaining part, which satisfies F2(q,r) = 0(fie 2 ). Recall that the functions 
afc(r) have zero average and therefore we have that 

(Fi> = 0. 

4.3 Existence of parameterizations of the perturbed invariant manifolds 

The next step is to prove the existence of parameterizations of the unstable and stable invariant 
manifolds of the periodic orbit given in Theorem 14.11 

As is done in [BFGSll] (see also |Bal06j ). we follow [LMS03| SauOlj . and we write the 
invariant manifolds as graphs of suitable generating functions which are solutions of a Hamilton- 
Jacobi equation in appropriate variables. To this end, we consider the symplectic change of 
variables (see |Bal06HBFG5II] ) 

f q = q {u) 

{p—T-y m 

I Po(u) 

where (qo(u),po(u)) is the parameterization of the homoclinic orbit given in Hypothesis HP6. 
This change is well defined for any u € C such that po(u) / and leads to a new Hamiltonian 
given by 

eH(u,w,t) = eH ( q (u), — ] , (32) 

where H is the Hamiltonian defined in ()29|) . When /x = 0, H becomes Hq defined in (|12j) . 
Then, the separatrix of the unperturbed system (/x = 0) for H can be parameterized as a graph 
as w = po(u) 2 . If we want to obtain parameterizations of the perturbed invariant manifolds, 
we can take into account the well known fact that, locally, they are Lagrangian and can be 
obtained as graphs of generating functions which are solutions of the Hamilton-Jacobi equation 
associated to the Hamiltonian eH. That is, we parameterize the invariant manifolds as graphs 
ast« = d u T u ' s (u,T), where the functions T u ' s satisfy 

d T T(u, t) + eH(u, d u T(u, r) , r) = 0. (33) 

The solutions of this equation give parameterizations of the invariant manifolds, which, in the 
original variables, read 

(q,P) = ?oW, r~, • (34) 



Po(u) 

The inconvenient of these parameterizations is that they are not defined everywhere since po (u) 
might vanish in some u & C. In some concrete happens for the classical pendulum, 

Po(u) 7^ for any u G C, and therefore we do not have any problem. However, this is not always 
de case. In particular, for the figure eight separatrices (see Figure [2]) there is always a real value 
of u such that po{u) = 0. For the graph case (see Figure [1]), po(u) ^ for all u £ M but it 
can vanish for complex values of the variable. Therefore, as was done in [BFGSll] . we look for 
parameterizations of the form (|34j) in very special complex domains for the variable u, which 
were defined in that paper and were called boomerang domains (see Figure ED . The reason of 
the somewhat strange shape of these domains is explained in full detail in that paper. Let us 
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just say here that we need domains which contain both an interval of the real line of width 
independent of e and points at a distance 0(e) of the singularities u = ±ia. Moreover, we want 
these properties to be still satisfied when we intersect the boomerang domains associated to the 
stable and unstable manifolds. These domains are defined as follows 

D S K d = {u £ C; |Imii| < tan (3\R,eu + a — ke, |Im u\ < tan ^Reu + a — ke, 

|Imu| > tan/^Reii + a — d} 
D" rf ={u£C; |Im-u| < — tan fiiHeu + a — ke, |Im u\ < tan/^Re u + a — ke, 

(35) 

|Imii| > tan ,$2 Re u + a — d} 

U {u £ C; |Imit| < — tan /?iReu + a — ke, |Imu| > — tan^Reu + a — d, 
Re u < 0} , 

where f3\ G (0, vr/2) is any fixed angle. The choice of /?2 goes as follows. First let us point out 
that the zeros of po(u) are isolated in C. Moreover, close to the singularities u = ±ia, po(u) can 
not vanish. Then, in order to assure that po(u) does not vanish in the whole domains D S K d and 
D™ d , one has to choose an angle such that $2 > Pi has a positive lower bound independent of 
e and /i and such that the lines |Imn| = tan ^Reu + a do not contain any zero of po(u). Then, 
taking e > and d > independent of e, both small enough, one can guarantee that po(u) does 
not vanish neither in D S K d nor in D^ d . 




Figure 3: The boomerang domains D™, and D S K . defined in ()35j) . 

As it is well known, to obtain the parameterizations of the invariant manifolds using equation 
(|33p . one has to impose the asymptotic conditions 

lim Pq 1 (u) ■ d u T u (u, r) = (for the unstable manifold) (36) 

Re u — > — oo 

lim Pq 1 (u) ■ d u T s (u,r) = (for the stable manifold). (37) 

Re u— >+oo 

These asymptotic conditions ensure that the invariant manifolds tend to the periodic orbit 
{q,p) = (0] 0) °f system associated to Hamiltonian ([29]) as Ren — > ±oo. Of course these 
conditions do not have any meaning in the domains D™' d since these domains are bounded. This 
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implies that, to prove the existence of the parameterizations of the invariant manifolds in these 
domains, one has to start with different domains were these asymptotic conditions make sense 
and then one has to find a way to extend them analytically to the domains D taking into 
account that these parameterizations become undefined at the points where po(u) = 0. This 
process is explained in great detail in [BFGS11 . Following the notation of that paper, the case 
we are considering now ressembles the case r = 1 and £ = 2 of that paper. Therefore here we 
just sketch it and give the main ideas. 

The main steps of this procedure are the following. 

1 We first obtain parameterizations of the form f)34[) satisfying the asymptotic conditions 
([36]) and (|37|) in the domains 



D^ p = {ueC;Reu<-p} 
D S ao,p = {u€C;Reu> p}. 

Note that thanks to the exponential decay of po(u) as Ren — > ±oo, one can ensure that 
for p big enough, po(u) ^ for u £ D^ p . 

2 To reach the boomerangs domains D™' d , we cannot extend the just obtained parameteriza- 
tions since in the analytical extension procedure, we might reach points where poiu) = 0. 
Thus, we switch to a different kind of parameterizations of the invariant manifolds, which 
are chosen of the form 

(q,p) = (Q(v,r),P(v,T)) (38) 

in such a way that (Q(v + es,t + s),P(v + es,t + s)) are solutions of the differential 
equation associated to Hamiltonian (|29p . We will extend these parameterizations and once 
we reach the domains D^' d , in which po(u) ^ 0, we will switch back to parameterizations 
of the form (|34p . This means that we need to define the parameterizations (|38p in domains 
which overlap both with the domains D^ p and with the domains D ' .. We define the 
following domains (see also Figured]), 



p' 



" = | u £ C; |Imii| < — tan /3iReii + a — ke, Re u > —p, 
|Imii| < — tan ^Reu + a — — | 
0°*'^ = jii G C; |Im u\ < tan /?iRe u + a — ke, Re u < p 
\lmu\ < tan /^Reu + a — - j, 



(39) 



where p' > p so that D° u d '* n Dlo )P / 0, * = u, s. 

The parameterizations fl38|) were used in [DS921 IDS971 IGel97al IGeWOI [BFOil IBFOK] . It is 
easy to see |Gel97a| that (Q, P) have to satisfy 

Q\_( P + fidpH^Q,^) 



Ce \ p J V -(V'(Q + ^(r))-^(T)))-M(ft^T) J ' 
where C £ is the operator 

C-e = E~ 1 d T + d v 
and H\ is the Hamiltonian defined in (|30p . 
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Figure 4: The domains D° n d '™ and D° u d '^ defined in ([39)) . 

3 To make this analytic extension procedure with parameterizations of the type (|38p . we 
first derive parameterizations (Q u ' s (v , r), P u ' s (v, r)) in the domains n D%o tP and 

D™^' s n p from the parameterizations (|34p obtained in Step 1. Taking into account 
the change of variables (|3ip . it is natural to look for these parameterizations as 

Q u > s (v,t) =q (v+U u > s (v,T)) 

= d u T^(v+U u > s (v,T)) 
l ' ' Po (v+U^(v,t)) ' 

where U u ' s define changes of variables u = v +U u ' s (v,t) so that (Q u,s , P u > s ) satisfy the 
system of equations ([30]) . 



4 Once we have obtained these parameterizations, we use equation (|4Up to extend them to 
the whole domains D ^'™ and D° p n ^- This step is straightforward since these domains are 
not close to the singularities u = ±zct. 

5 Now we have parameterizations of the invariant manifolds of the form (|38p in the following 
domains, which are depicted in Figure El 

1 K,d ~ U p,d,K 11 U K,d . . 

,-out,s _ n out,s _ nS v / 

1 K,d - u p a^ 1 1 M' 

where, by construction, po(u) does not vanish. Thus, we can use these domains as transi- 
tion domains where we can go back to the parameterizations (|34p and where the Hamilton- 
Jacobi equation ([33]) can be used. To obtain them, we look for changes of variables 
v = u + V u ' s (u, r) which satisfy 

Q^( u + V u > s (u,T),T)=q (u), 

where Q u ' s are the first components of the parameterizations just obtained. From these 
changes, we can deduce the generating functions T u ' s which give the parameterizations 
(GT 
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Figure 5: The domains I° u d ' u and I° u d ' s denned in (HU- 

6 The last step is to extend analytically the parameterizations (|34p just obtained to the whole 
domains D™'*,. Since in these domains po(u) 7^ 0, we take into account that the associated 
generating functions T u,s satisfy equation (j33j) and we use this equation to obtain these 
final analytic extensions. 

Note that in the case that po{u) 7^ for any u G C, as happens for the classical pendulum, one 
can skip steps 2 to 5. That is, one can extend the parameterizations (|34p from infinity directly 
to the boomerang domains D™'% 

As a conclusion of this procedure, which is explained in fully detail in Sections 6 and 7 of 
|BFGSll] . we can state the following theorem. Recall that using the notation of that paper, the 
problem considered in the present paper is stronly related to the case r = 1 and £ = 2 considered 
in that paper. We state the result for the unstable invariant manifold. The stable one satisfies 
analogous properties. 



Theorem 4.3. Let fj>o be the constant given in Theorem \4-l\ d\ > 0, k± > big enough and 
So > small enough. Then, for [i G B(fJ>o) an d £ G (0,£o), the unstable invariant manifold of 
the periodic orbit (q,p) = of the system associated to Hamiltonian (|29p has a parameterization 
of the form (|MD for (u,t) G D u Ki di x T CT . 

Moreover, there exists a real constant b\ > independent of e and [i such that, for (u, r) G 
di x T CT , the generating function T u which gives the parameterization ()34|) satisfies that 

\d u T u (u,r)-d u T (u)\< 

\u z + or 1 



where 



a 



T (u)= / p 2 (v)dv, (42) 
J — 00 

is the parameterization of the unperturbed separatrix. 

4.4 The asymptotic first order of d u T u,s close to the singularities ±ia 

Theorem 14.31 gives parameterizations of the invariant manifolds up to points at a distance of 
order 0{e) of the singularities. Nevertheless, at this distance of the singularities Theorem 14.31 is 
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not a perturbative result with respect to the singular parameter e. Namely, at a distance 0(e) of 
u = Ha the parameterizations of the perturbed invariant manifolds are not well approximated 
by the parameterization of the unperturbed separatrix. This fact can be easily seen from (|15p 
and the bounds given in Theorem 14.31 since for u =F ia ~ e both 8 u Tq(u) = Pq(u) and the 
remainder become of size 0(l/e 2 ). Then, to study the difference between the manifolds, we 
need to look for better approximations of T u,s in the following domains, which are usually called 
inner domains (see Figure [6|), 

^ >1 k,c~' U = {it G C; Imu > — tan f3\ (Re u + ce 1 ) + a, Im u < — tan /^Re u + a — ke, 

Im u < — tan /3oRe u + a — ke} 
D^ = {uGC;uGD^} (43) 

D™y s ={u£C;-ueD^ u } 

where k > 0, c > and 7 G (0, 1). The angles ft\ and 02 are the ones considered in the definition 
of the boomerang domains in (|35p and (3q is any angle satisfying that fi\ — /3q has a positive 
lower bound independent of e and /i. Let us observe that, if u £ D 1 ^ '*, * = n, s, and e is small 
enough, then O(ke) < \u =F < 0(e 7 ) and u G D* d , * = u, s. 




— ia 



Figure 6: The inner domains defined in (|43p. 

We obtain good approximations of the functions T u,s in these domains through a singular 
limit. Since the study of both invariant manifolds close either to u = ia or u = —ia is analogous, 
we only study them in the domain D^' c Thus, we consider the change of variables 

z = e _1 (u — ia). 

The variable z is called the inner variable, in contraposition to the outer variable u. By definition 
of To in (|42p and using the expansion around the singularities of po(u) in (|15p . we have that 

<9 u T (ez + ia) = p 2 {ez + ia) = ^(l + ((ez) 2 / M )) 
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and, using the results of Theorem 14, 3\ we have that 



\d u T u > s {ez + ia,r) - d u T (ez + ia)\<K 



e 2 |z| 3 



Hence, to study the first order in e of the functions T u,s at a distance 0{e) of the singularities 
we scale the generating function as 

^ s (z, t) = eC+ 2 T u ' s (ia + ez, r). (44) 

Then, the Hamilton-Jacobi equation (|33p becomes 

d T ?p + e 2 C+ 2 H (ia + ez, e~ 2 C\d z ^, r) = 0, (45) 

where H is the Hamiltonian function defined in (|32|) . The corresponding Hamiltonian is 

H(z, w, r) = e 2 C+ 2 H (ia + ez, e~ 2 C\w, r) . (46) 

We study equation (|45p in the domain Pj^c ' u x T CT , where 

={z€C;za + ez€ D^} . (47) 



To study equation (|45p . as a first step it is natural to study it in the limit case e = 0, which 
reads 



3 T ^o + \z 2 (d z iPof 



\ Z Z I 



where 

-2 



G{w) = C+ A G(C+w) 

The other terms are the ones coming from the functions F and R defined in ([9]) and (|10p 
respectively. Indeed, fm(r)/z 2 is the leading term coming from the degree M terms of F(x,t) 
(recall that F is a trigonometric polynomial in x of degree M thanks to Hypothesis HP2), and 



Z t>l 

comes from the degree M terms of the function R(x, I, r) in (jlOD , 

Equation (|48p is the Hamilton-Jacobi equation associated to the non-autonomous Hamilto- 
nian 

1 



~H (z,w,t) = t;Z 2 w 2 



^2 + » ^(zw) + + ^2 E ^ j • 



which satisfies that T~L — > T~Lq as e — > 0, where % is the Hamiltonian function defined in (|46p . 

Remark 4.4. ^4s an example we show the inner equation associated to the Hamiltonian (|2l 
considered in Section {K4\ Recall that, as we have seen in Section \2Ji\ the integrable system Hq 
associated to (f24"j) is simply the classical pendulum, whose upper separatrix has singularities at 
u = ±m /2km , k £ Z and C+ = — 2z. Then, the inner equation around u = in/2 is given by 

dzi>o + \z 2 (d z ^ ) 2 ~ 2^(1 + sinr) M - ^ ) ~ 2 V i (zd^of = 0. 
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Figure 7: The domains V K 'g and V K g defined in (@9j) . 

In |Bal06j , the author studies an inner equation slightly different from fl35J) (see also |OSS03j ) . 
The main differences between the inner equation considered in |Bal06] and (|48|) is that in |Bal06| . 
the author assumes G = and assumes also polynomial dependence on d z i/jo whereas fi8j) can 
have the full jet. Nevertheless, one can easily see that using the techniques considered in |Bal06j 
one can obtain exactly the same results obtained in that paper for inner equations of the form 
(|48|) since G is cubic and therefore, for z big enough, is smaller than the other terms. We state 
these results in the next theorem, but first we introduce certain domains. The solutions of (|48p 
are studied in the complex domains 

V+£ = { Z £C; |Im z\ > 9 Re z + k} 

« = Hc;-z £ d;;} (49) 

for k > and 9 > 0. Let us observe that, for any c > 0, U^c '* C D^* n ^ 2 for * = u,s. 

Nevertheless, since through the proof we will have to change the slope of the domains "D C, we 
start with a slope 9 = tan/32/2. The difference between the stable and unstable manifolds of 
the inner equation is studied in the intersection domain 

Kb = V+ J n V tfi n i z e C; Im z < 0} . (50) 



Theorem 4.5. For any fx G B(^,q) the following statements are satisfied: 

1. There exists K2 > such that equation (|48|) has solutions tpQ : ^^\ an p 2 /2 x T ff -> C, 
* = u, s, of the form 



^ s (z,r) = — + ^ s ( z ,t)+K u > s , K u > s G C 



(51) 



where tp ' are analytic functions in all their variables. Moreover, the derivatives of ip ! 
are uniquely determined by the condition 



sup 



2 , tan /3 2 /2 



z 3 d z %(z,T) 



< oo 
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—in 




Figure 8: The domain H^o defined in ([50|) , 



for * = u,s. In fact, one can choose ip ' such that 



sup 



2, tan /3 2 /2 



xT 



< oo 



for * = u,s. 

2. There exists K3 > K2, analytic functions {x^(/ u )}/ og z- defined on B(hq) and ( : TZ^ s tan ^ 2 x 
T a — > C such that two solutions iPq' s of equation (|48p 0/ the form given in (|5ip u>ii/i 
K u = K S , satisfy 

- V>g) (*,r) = ^^x [fc] (/^)e lfc(2 - T+K(2 ' r)) . 

fc<0 

Moreover, the function £ is 0/ t/ie /orro 

C(^,r) = -nblnz + ((z,t), 
where b is the constant defined in (|17p and £ satisfies 



sup 

v ' ' Kg ,tan P2 



IzC(z,t)| < °°- 



The proof of this theorem follows the same lines as the theorem proved in |Bal06| . Moreover, 
it can be easily seen that the analytic functions {x^ (v)} k^z- are entire. 

To have a better knowledge of the parameterizations of the invariant manifolds in the inner 
domains X>™c '*> * = u, s in (|47p . we need to compare the parameterizations ip u,s , which are 
solutions of (|45p with iPq' s which are solutions of (|48p and have been given in Theorem 14.51 
Recall that Z>™c +1 * C ^^tanfe' * = u > s > an d therefore V>o' s are defined in these domains. 

We state the next theorem for the unstable invariant manifold. The stable manifold satisfies 
analogous properties. 
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Theorem 4.6. Let 7 G (0, 1), the constants K\ and K3 defined in Theorems \4-3\ and \4.5\ c% > 
and £q > small enough and K4 > max{/ti, K3} 6i<? enough, which might depend on the previous 
constants. Then, for e G (0,£o) an d A* G B(hq), there exists a constant hi > smc/i i/iat /or 

(^r)ePi n '+; u xT a) 



\d z r(z,r)-dMz,r)\ < 



b2£ M 

k M 



where 7 enters in the definition ofD^d™, V'o ^ s given in Theorem \4-5{ ip u is the scaling of the 
generating function T u given in (|44p and M /ias 6een given in HP2. 



The proof of this theorem follows the same lines as the proof of Theorem 4.16 of [BFGS11]. 
4.5 Study of the difference between the invariant manifolds 



Once we have obtained parameterizations of the invariant manifolds of the form (|34|) in the 
domains D S K . and D" . and studied their first order approximation near the singularities, 
the next step is to study their difference. To this end, we define 



A(u,r) = T u (u,t) -T s (u,t) 

in the domain R K 

Rn,d = {u £ C; |Im u\ < tan /^Re u + a — K£, |Im u\ > tan ,$2 Re u + a — d, 



D S K d R .D" d which is defined as 



|Im u\ < — tan /3iRe u + a — ke} 



(52) 



and can be seen in Figure 

We recall that po(u) 7^ if u G an d hence in this domain we can use the functions T s ' u . 



""--^ ia 


i(a — ke) 






/ --^K,d 










— ia 





Figure 9: The domain R K ( i defined in ([52]) . 

Subtracting equation (|33[) for T" and T s , one can easily see that A satisfies the partial 
differential equation 

C £ i = 0, (53) 
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where 

C £ = e~ 1 d T + {l + A{u,T))d u 

with 

1 

2pf(u) 



A ^ T ) =1727^ (duT^(u,T) + d u T{(u,r)) 



+ — rr / «piJi q (u),po{u) H — ,r ds, 

where #1 is the function defined in dH and T u ' s (u, r) = r (u)+2^ ,a (it, r) with <9 u T (n) = 
and T"' s are given in Theorem 14.31 

Following [Bal06j (see also |BFGSlf] ). to obtain the asymptotic expression of the difference 
A, we take advantage from the fact that it is a solution of the homogeneous linear partial 
differential equation ([53]) . In |Bal06] it is seen that if ([53]) has a solution £o such that (£q(u, t), r) 
is injective in R K d x T CT , then any solution of equation (|53p defined in R K) d x can be written 
as £ = T o £g for some function T. 

Following this approach, we begin by looking for a solution of the form 

£o{ u i T ) = e ~ u — t + C(u, t) 

where C is a function 27r-periodic in r, such that (£o( u > f), r) is injective in R K ,d x T CT . Moreover, 
to study the look for the first order of A, we need to compare such function C with the function 
£ obtained in Theorem 14.51 in the inner domains (|43|) . 

Theorem 4.7. Let us consider the constants d\ > defined in Theorem \4-3\ and K4 > in 
Theorem \4-6[ Sq > small enough and K5 > K4 enough, which might depend on the previous 
constants. Then, for e G (0, £0), /i G B{[Xq) and any k > K5 snc/i i/iai £k < a, i/iere exists a 
real-analytic function C(u,r) : R K di x T CT — > C snc/i i/tai 

• £o(w, t) = £ _1 n — r + C(u, r) is solution of ([53]) and 

(Co(«, r),r) = (e _1 ii - r + C(n, r), r) 

is injective. 

• There exists a constant 63 > independent of n, e and k, such that for (u, r) G R K ,d\ x ^a, 

\C(u, t)\ < 63 |//| In \u 2 + a 2 | 
|a»C(u,r)|< ..^W - 



|ir + 

Moreover, if we consider the constant c\ given in Theorem |^.6] and 7 satisfying 

M 



< 7 < 1, 



2 + M 

where M has been defined in Hypothesis HP2. Then, there exist a constant C((jl) defined 
for [i G B(fio) depending real- analytically in [i and a constant 64 > such that \C(p)\ < 

bi\n\ and, if (n,r) G (D 1 ^' u n £>™c|' s ) X T a for any k > k 5 , 



\C(u, t) - C(p) + fiblae - (e _1 (n - ia), r) | < 



it — ia 
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where £ is the function given in Theorem \4-5\ and b is the constant defined in (fTT|) respec- 
tively. 

The proof of this theorem follows the same lines as the proofs of Theorem 4.21 and Propo- 
sition 4.22 of [BFGSllj . 

As we have explained, since A = T u — T s is a solution of the same homogeneous partial 
differential equation as £o given in Theorem 14. 71 there exists a function T such that A = To^, 
which gives 

A(u, r) = T (e^u - r + C(tt, r)) . 

Since A is 27r-periodic in r, the function T is also 2-7r-periodic in its variable. Therefore, consid- 
ering the Fourier series of T we obtain 

A(n, t) = Y^ T^e ik ^ lu - T+C{u ' r) ) . (54) 
fcez 

Now we find the first asymptotic term of A which is strongly related with (tpQ — i/jq) (e" 1 (u — 
ia),r), where iPq' s are the solutions of the inner equation given in Theorem 14.51 We introduce 
the auxiliary function 

A+(u,r) = ^r [ * ] e ik ( £ - lu - T+c ^ T) ) 

k<0 

with 

T [ fc ] = ^ x ik]^ )e - lJ P-im-c^+^e) (55) 

where {x k (f J -)} k<0 are the coefficients given in Theorem 14.51 and C(^) and b are the constants 
obtained in Theorem 14. 71 and Theorem 14.51 respectively. The scaling C\je comes from the inner 
change in (|44p . 

We also introduce 

A- («,r) = Y^rfe ik ( £ - lu - T+c ^) 

k>0 

with 

e 

The function Aq (u,t) corresponds to the difference of the solutions of the inner equation close 
to u = —ia. Taking t,/i£1, Aq is nothing but the complex conjugate of Aq . In fact, as we 
know that A is a real analytic function in the u variable for real values of /i, r, we can define 
Aq as the function that satisfies that Ao = Aq" + Aq is also a real analytic function t,/i6R. 
We will see that the first order of A is given by 

A (u, r) = A(j"(n, r) + Aq (it, r), 

which can be written as 

A (u,r)= Yl T lk] e lk ( £ ' lu - T+C ^), 

k£Z\{0} 

where Y ^ are defined either by (|55[) and (|56p . For convenience we introduce T °' = 0. From 
now on, we consider real values of r £ T ff nl. 



T [k] = WM-[- fe ] (/i)e -^ +i | fc |(-C(M+M61n £ )_ (56) 
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Theorem 4.8. Let us consider the mean value ofT, Y^, defined in (j54p . s < 2/M, and Eq > 
small enough. Then, there exists a constant 65 > such that for e G (0,£o) and (i G B(fJ-o) CI M 
and (it, r) G (-R s in(i/e),di Hi) x T, £/ie following statements are satisfied. 

A(u,t)-T' ' - A (u,r 
|<9 u A(u,t) -0 M A o (u,r)| < 
|^A(«,r)-^Ao(«,r)| < 



< 


h\tA - 






e\ lne| 




< 




-^+/ilm 6 In e 




e 2 lne 




< 


£»5 ImI 


- ^ +/zlm 6 In e 




e 3 lne 





The proof of this Theorem follows the same lines as the proof of Theorem 4.23 in [BFGSllJ. 
Note that, following the notation of (BFGSllJ, the case considered in Theorem 14.81 corresponds 
to r = 1 and £ = 2 in |BFGSllj and thus £-2r = 0. 

We observe that 8 u Aq gives the correct asymptotic prediction of d u A if Yq ^ 7^ 0. In fact, 
we only need this coefficient to give a simpler leading term of the asymptotic formula. For this 
purpose let us define the function 

f^) = Clx [ - 1] (At), 

where C+ is the constant defined in (|15p and ^(m) ^ s the constant given in Theorem l4.51 This 
function /(//) is the one appearing in the asymptotic formula for the area of the lobes (|16p given 
in Theorem 12.51 From its definition just given, we can see that this constant /(/x) is essentially 
given by x}-~ 1 \ij), that is, by the first harmonic Fourier coefficient of the difference between the 
solutions iI>q' s of the inner equation ()48|) (see Theorem I4.5() . In particular, the zeros of f(fJ-) 
correspond to the zeros of x'~^(aO- 
We define 

A 00 (n,r) = ^e^Re (/(^e-^"— ^^(f-^r))^ 

where the constant b has been defined in (|17p and the function C and the constant C(/i) have 
been defined in Theorem 14.71 

Corollary 4.9. There exists a constant b§ > such that for e G (0,eo), A* G B(/j,q) n K and 
(it, r) G (-ff s ln(l/e),<ii H R) x T, the following statements are satisfied, 



A(u,t)-T^-A 00 (u,t) 



&6 \lA --+^Imfeln £ 



e lne 



|d n A(n,r)-d M A o(u,T)| < _^L e -f+Mim6lu e 

e z | lne I 

kj2 A / \ o2 a / \\ ^ ^6 H --+ulm61ne 
\ d u A { U ^ T ) ~ d u A 00 (u,T)\ < ^r- re ^ 

1 e d |me| 

To finish the proof of Theorem 12.51 it is enough to proceed as in Section 4.8 of [BFGSli] . 
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